Asteroids, Comets ; Meteors 1991, pp. 605-608 
Lunar and Planetary Institute, Houston, 1992 


■5VV/-/3 


605 



DETERMINATION OF ORBITS OF COMETS: P/KEARNS- 
KWEE, P/GUNN, INCLUDING NONGRAVITATION AL 
EFFECTS IN THE COMETS’ MOTION 

B. Todorovic-Juchniewicz, G. Sitarski 

Space Research Centre, Bartycka 18, 00-716 Warszawa, Poland 

ABSTRACT 

To improve the orbits we collected all the positional observations of the comets. The 
observations were selected and weighted according to objective mathematical criteria and 
the mean residuals a priori were calculated for both comets. We took into account non- 
gravitational effects in the comets’ motion using the Marsden’s method applied in two 
ways: either determining the three constant parameters A\, A 2 , A 3 or the four parameters 
A,rj,I,<f> connected with the rotating nucleus of the comet. To link successfully all the 
observations we had to assume for both comets that A(t) = A$ exp (— B ■ t) where B was 
an additional nongravitational parameter. 

METHOD FOR NONGRAVITATIONAL EFFECTS 

According to the known Marsden’s method a nogravitational force acting on a comet 
is defined by its three components in orbital coordinates: F\ - the radial component, Fi - 
the transverse component in the orbit plane, F 3 - the component normal to the orbit plane 
(Marsden 1969). It is assumed that Fi = A,g(r), and the nongravitational parameters A,, 
i = 1,2,3 are to be determined as constant quantities from observations for an individual 
comet; the analytical form of the function g{r ) is known (Marsden et al. 1973). 

If the nongravitational force acts on the comet owing to the rocket effect of ejection 
of material from the rotating comet’s nucleus, then: 

Ai = ACi(r},I, A = v + <f>) 

where A = ( A\ + A\ + Af) 1 / 2 , 77 is the lag angle, I - the equatorial obliquity, v - the true 
anomaly of the comet, <j> - the cometocentric solar longitude at perihelion; the direction 
cosines Ci have the following form: 

Ci = cos 77 + ( 1 — cos 77) sin 2 I sin 2 A 

C 2 = sin 77 cos I + ( 1 — cos 77 ) sin 2 J sin A cos A 

C3 = — [ sin 77 cos A — ( 1 — cos rj ) cos I sin A] sin I 

It appears that A , 77, 1, <f> can be determined from observations by the least squares 
method (Sitarski 1990). To link successfully all the observations of P/Kearns- Kwee and 
of P/Gunn as well, we had to assume that the nongravitational parameter A depended 
on time. We accepted that A(t) = Ao exp (-B • t ) where t was counted in days from the 
epoch of osculation. Thus the five nongravitational parameters: Ao,B,r),I,<j) have been 
determined from observations along with corrections to the six orbital elements. 
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OBSERVATIONAL MATERIAL 


Observations of the comets were s electe d and weig hte d (each apparition separately) 
according to objective mathematical criteria (Bielicki 1972). A mean residual o priori was 
calculated as follows: mean residual = [£*«., u> ; 7 z^/E*nj] 1/2 where ^-number of appari- 
tions, n ; - number of residuals of the j-th apparition, wj — weight oTthe j'-th apparition, 
fij — value of the mean residual of the j-th apparition. •• 

Characteristics of observations of Comet P/Kearns-Kwee: 


Apparition 

Number of 
observations 

Number of 
residuals 

Mean 

residual 

Weight of 
apparition 

1963.08.17 - 1965.04.24 

60 

115 

1"33 

0.64 

1971.07.26 - 1973.04.04 

96 

188 

1.26 

0.71 

1981.06.29 - 1982.02.26 

46 

83 

1.41 

0.56 

1989.09.10 - 1991.03.21 

64 

123 

0.74 

2.07 


Observation interval: 1963 Aug. 17 — 1991 Mar. 21 

Number of observations: 266 

Number of residuals used for the orbit improvement: 509 

Mean residual a priori : l!'06 

In case of Comet P /Kearns Kwee we detected a displacement of the photometric cen- 
ter from t he cent er of mass of the comet’s nucleus (Sitarski 1984). We assume that the dis- 
placement of both centers along the radius- vector r is expressed by Ar = Dr' 3 exp(-r 2 /2) 
where the value of D is to be determined from the observational equations by the least 
squares method along with other parameters of the comet’s motion. 

Characteristics of observations of Comet P/Gunn: 


Interval 

Number of 

Number of 

Mean 

Weight of 


observations 

residuals 

residual 

interval 

1954.08.08 - 1954.08.08 

2 

4 

“ 

1.00 

1970.10.27 - 1971.11.18 

13 

26 

0"74 

0.98 

1972.12.03 - 1975.07.03 

42 

83 

1.07 

0.46 

1976.04.24 - 1978.10.28 

24 

46 

1.04 

0.49 

1980.12.07 - 1983.09.28 

20 

40 

0.81 

0.81 

1984.08.24 - 1987.01.04 

20 

38 

0.39 

3.50 

1988.01.15 - 1990.09.24 

78 

153 

0.78 

0.88 


„ r—. -r; _r- ■ : ; 

... . : 7 ; ^ - 


... ... 


Observation interval: 1970 Oct. 27 - 1990 Sep. 24 


Number of observations: 197 

Number of residuals used for the orbit improvement: 386 

Mean residual a priori : 0”73 
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RESULTS 


Orbit of Comet P/Kearns-Kwee: 

(Epoch of osculation: 1963 Nov. 27.0 ET, Equinox: 1950.0) 

1. With nongravitational parameters Ai,A 2 , A 3 : 

T = 1963 Dec. 7.00738 ET u> = 131°19146 
q = 2.21317774 a.u. Cl = 315.43935 

e = 0.48653298 *= 8.99170 

A ! = (+0.59457 ± 0.13152) x 10~ 8 
A 2 = (-0.38585 ± 0.00154) x 10~ 8 
A 3 = (-0.15898 ± 0.20130) x 10" 8 
D = (+0.95212 ± 0.17450) x 10~ 2 

Mean residual a posteriori : 4”67 

2. With A = Aoe~ B t and angular parameters r 7, I , <f> : 

T = 1963 Dec. 7.00735 ET u = 131°19250 

q = 2.21318442 a.u. Cl = 315.43918 

e = 0.48653473 * = 8.99151 

+0 = (+1.47950 ± 0.06238) x 10~ 8 
B = (+0.11971 ± 0.00151) x 10 -3 /day 
77= 33°30± 2° 68 

I = 143.22 ± 8.55 

(f> = 334.45 ± 17.01 
D = (+0.36178 ± 0.04559) x 10~ 2 

Mean residual a posteriori : 1“ 24 


Before the discovery of P/Kearns-Kwee a close aprroach of the comet to Jupiter to 
within 0.033 a.u. occured in November 1961. That approach considerably changed the 
comet’s orbit (before the approach: q = 4.302, e = 0.691, i = 2°88). 

Orbit of Comet P/Gunn: 

(Epoch of osculation: 1989 Oct. 1.0 ET, Equinox: 1950.0) 

1. With nongravitational parameters +1 , +.2 , ^3 : 

T = 1989 Sep. 24.96595 ET oj = 196°93927 
q = 2.47155291 a.u. Cl = 67.86640 

e = 0.31439342 »' = 10.37272 

Ai = (+2.52440 ± 0.00545) x 10" 8 
A 2 = (+0.50284 ± 0.00050) x 10" 8 
A 3 = (-0.14875 ± 0.15015) x 10“ 8 
Mean residual a posteriori : 1"52 
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2. With A = A 0 e~ B t and the angular parameters : 

T = 1989 Sep. 24.97004 ET w = 196*94062 
q = 2.47155081 a.u. D = 67.86623 

e = 0.31439419 i = 10.37270 

A 0 = (+1.15530 ± 0.07298) x 10" 8 
B = (+0.12002 ± 0.00751) x 10- 3 /day 
rj — 53*45 ± 12*91 
I = 73.29 ± 3.67 

4> = 84.84 ± 4.37 

Mean residual a posteriori: lf ; 15 

Comet P/Gunn approached Jupiter to within 0.353 a.u. in July 1965, and it changed 
the comet’s orbit (before the approach: q = 3.290, e = 0.178, i — 10. 87). 


CONCLUSIONS - 

We can see that accepting the constant nongravitational parameters Ai,A 2 ,A 3 it 
was impossible to find resonable solutions for orbits of the considered comets, the mean I 

resid ual s a posteriori were much greater than those a priori (especially for P/Kearns- j 

Kwee). Assuming an expotential time dependence for the nongravitational parameter 
A(t) we got the acceptable solutions for both comets. 

Both comets approached Jupiter, and the approaches changed the comets orbits in 
the same manner: perihelion distances diminished. The expotential decreasing of non- 
gravitational effects, as detected in the motion of both comets (the positive values of the 
parameter B ), can be explained by the diminishing of the comets’ activity with time after 
a rapid increase of activity caused by the changes of orbits due to approaches of the comets 
to Jupiter. 
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